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We investigate the statistical mechanics of glass-forming materials and jamming matter by means 
of a geometrically driven approach based on a revised cell theory. By considering the system 
as constituted of jammed blocks of increasing sizes, we obtain a unified picture that describes 
accurately the whole process from low densities to limit densities at the glass/jamming transition. 

The approach retrieves many of the aspects of existing theories unifying them into a coherent 
framework. In particular, at low densities we find a free volume regime, based on local relaxation 
process, at intermediate densities a cooperative length sets in, where both local and cooperative 
relaxation process are present. At even higher densities the increasing cooperative length suppresses 
the local relaxation and only the cooperative relaxation survives characterized by the divergence of 
the cooperative length, as suggested by the random first order theory. Finally a relation between 
the cooperative length and the hyperuniform length is also suggested. 


There are many systems, such as molecular liquids, colloids, granular materials, foams and others that by changing 
the control parameters exhibit a very slow dynamics followed by a structural arrest or quasi structural arrest. Although 
many progresses have been made, understanding the glass transition, or jamming, from a fluid phase to an amorphous 
solid is still one of the major problems in condensed matter. In particular, the quest concerns understanding the 
main mechanism for the huge increase of the viscosity, or relaxation time, as the temperature decreases or the density 
increases that is observed across different systems. 

Many theories have been proposed. Among the first ones is the Free volume theory mostly developed by Cohen 
and Turnbul m, and the Cooperative rearranging regions (CRR) of Adam and Gibbs [1]. Both, using simple and 
physical intuitive arguments, predict a strong exponential divergence at the ideal glass transition characterized by a 
Voghel-Fulcher-Tamman (VFT) law. Free volume theory is based on a local property of the particles: each particle 
can move in a free volume where the center of a particle can translate, given that all other particles are fixed. Roughly 
we can say that the particle rattles in a cage of its free volume. At very high density, or at low temperature, close 
to the glass transition, relaxation occur when very rarely a particle ends up in a cage with large enough free volume 
and manage to jump out of the cage. A somehow opposite view was assumed by Adam and Gibbs, who argued that 
relaxation occurs due to a cooperative rearrangement of group of particles. The linear size of such region ^ diverges 
at a temperature Tg, the Kautzmann temperature (3, where the configurational entropy vanishes. 

An important advance in the theory of the glass transition is represented by the Mode coupling theory mostly 
developed by Goetze and collaborators [i[7]. Starting from first-principles they derived an equation for the time 
dependent density autocorrelation function, which makes a number of precise dynamical predictions. In particular, 
they predicted that the relaxation time should diverge as a power law at a critical temperature, where structural 
arrest is reached. However, molecular glasses do not show such power law divergence and, nowadays, Tc is considered 
a crossover towards a hopping regime characterized by exponential type of the relaxation time. The nature of the 
approximations involved in the mode coupling theory, has been elucidated by the spin glass theory which predicts, in 
mean held approximation, the same solution, in its schematic version, as the mode coupling theory. 

In recent years Random First Order Transition (RFOT), hrst introduced by Kirkpatrick, Thirumalai and Wolynes 
[5] and later developed by Wolynes and collaborators 0 [TU], has become one of the most popular theory of the glass 
transition, which, besides the idea of the cooperatively rearranging regions, combines also features of spin glass m 
and mode coupling theory. It predicts an apparent power law divergence of the relaxation time, followed by a Voghel- 
Fulcher-Tamman law diverging at the Kautzmann temperature. Each theory captures some correct aspect of the 
glass transition. It is then plausible that the ”hnal” theory must reproduce the essential features of these different 
approaches. 

In this paper, we show that a revised version of the Gell Theory of the Glass Transition [12] reproduces the essential 
ingredients of both Free Volume and Cooperative Rearranging Regions in a unified manner leading to new predictions. 

Free Volume Theory. Let N{n) be the number of particles with free volume Vf{n), where n is a discrete index. 
The free volume distribution can be calculated by maximizing, VL{N(n)), the number of ways of redistributing the 
free volume, Vl{N{n)) = V(n)! under the following constraints X)where 

Vf is the total free volume and N is the number of particles. 

After maximizing Q{N{n)) and passing to the continuum limit one obtains 
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where p(vf)dvf is the probability to find a particle in a free volume between Vf and Vf + dvf and (vf) = Vf/N is the 
average free volume. In this theory the relaxation time r is related to the inverse of the probability that a particle 
jumps out of the cage. If a particle is in a cage with small free volume Vf the probability to escape Pesc will be small. 
In a simple approach it is assumed that, below a threshold Vf < Vj, the particle is localized, otherwise it can jump 
and get out of the cage. In this case the probability Pesc that a given particle jumps out of the cage is given by the 
probability that the particle has a free volume Vf > Vj : 


Pesc = / p{vf) dvf , ( 2 ) 

Jv) 

where p{vf) is given by Eq0 The relaxation time t/tq (where tq is a microscopic time) is proportional to the inverse 
of this probability: 


— = = exp{v}/{vf}) , (3) 

'^0 ^esc 

Within the free volume theory, Cohen and Turnbull estimated that the average free volume goes to zero linearly with 
the temperature 


{vf)=AiT-To) , (4) 

Inserting Eq|^ into Eqj^ gives a VET law for the relaxation time 

- =exp(i?/(T-To)) , (5) 

To 

where A and B are constants. 

Although the relaxation time for many glass formers can be reasonably well fitted by a VET law, the prediction of 
a simple exponential for the free volume distribution, EqlTJ is not supported by numerical simulations on hard spheres 
and experimental data on granular materials [IlHIl] which instead show a very good fit with a Gamma distribution 
(see later Eqj^. 

Cooperative Rearranging Regions. An alternative theory was introduced by Adam and Gibbs in 1958 [1], who 
introduced the concept of cooperative rearranging regions. 

The main idea underneath this approach is that, close to the glass transition, due to the crowding of the particles, 
the decay towards equilibrium of a density fluctuation is due to a cooperative rearrangement of an entire region. A 
cooperative rearranging region of linear size ^ can be defined as the smallest region that can be rearranged without 
involving particles outside its boundary. It is argued that the relaxation time diverges exponentially with the excess 
entropy with respect to the underline crystalline state Sc- 

— =exp{C/{Tsc)) , (6) 

To 

where C is a constant. Adam and Gibbs argued that Sc is proportional to with d the space dimension. Following 
Kauzmann suggestion that the excess entropy goes to zero linearly at Tk, it follows that the relaxation time diverges 
according to the VET law Eq|3]with Tq = Tk- 

Wolynes and co-workers [^ . ^ fTO] generalized the Adam and Gibbs approach, proposing a mosaic picture, where 
the system is partitioned in droplets of linear size Like in ordinary first order transition the free energy of nucleating 
a droplet contains a volume term plus a surface term. Within the random first order transition the surface term is 
modified into a term which is proportional to with 9 = d/2, (see also [TB] and Berthier and Biroli m for further 
elaboration and a discussion of other possible values of 9). Gonsequently the configurational entropy results to be 
related to the size of the droplet as Sc = const/T//^~^ and the relaxation time is 


T 
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where AF is the droplet reconfiguration free energy barrier and it is given by 

, ( 8 ) 

Given that 9 = d/2 the Adam and Gibbs relation, Eq|6| is recovered. 

Cell Theory. In [TUlIHlIini the authors of the present paper, based on lattice theories of liquids[2n]) introduced a 
‘cell theory’ that combines the ideas of inherent structures, free-volume theory and geometrical packing properties to 
derive a general theory to understand the complex dynamics of glass-forming liquids, granular packings and amorphous 
solids. The main feature of this theory is the demonstration that thermodynamical properties of these systems can 
be retrieved from the study of geometrical and topological properties of local configurations only. 

In the cell theory, the partition function of a system of N particles in a volume V is: 

Z = . (9) 

{N{n)} 

where 

F({A(n)}) = ^ A(n){e(n) - fcr[ln^^^-lniP(n)]} , (10) 

n 

is the free energy [3T] which is dependent on the distribution of the cell-shapes and sizes {A(n)} [TH UHl HH] ; ^^/(n) 
is the ‘free volume’ associated with a particle in a cell with a set of geometrical and topological parameters n; e(n) 
is the energy associated with a particle in a cell with n; V{n) is associated with the probability to find a cell with 
n that is not single-occupied [50]. The quantity r2({A(n)}) in Eqj^ counts the number of distinct space-partitions 
(associated with the inherent states) made with the same set of {A^)}. 

The key elements to be estimated in Eq{T^ are n({A(n)}) and V{n). At high density particles move only locally 
and all Voronoi cells are singly occupied and 'P{n) ^ 1. For what concerns the terms n({A(n)}), it was pointed out 
in miiisKisi that the maximum number of distinct configurations that -in principle- can be made by positioning in 
different ways the N cells distributed in groups of A^(n) is A!/J([^ A(n)!. However, this maximal value cannot be in 
general achieved by H since some of these combinations of cells do not generate space-filling assemblies and others 
might be not associated with any inherent state. As a crude expression of this we can say that only a fraction of 
cells N/X can be considered in this exchange and, analogously, for each kind of cell only the fraction N{n)/X can be 
exchanged. The idea follows from the mosaic picture of the cooperative rearranging regions discussed before that the 
system can be partitioned in a mosaic each region of linear size ^ and in each region the effective number of degree of 
freedom is proportional to therefore the total number of degree of freedom N is reduced by a factor 

A = , (11) 

with tq a characteristic size such that pr^ = 1, with p the density (the volume fraction). Under this assumption we 
have: 


H({iV(n)}) = 


(N/xy. 


n„(lV(n)/A)! ■ 


From Eqsj^l^ and[T2l we obtain 


Z = exp| -/3^ A(n)[e(n) - A:T( 

{N(n)} n 


In 


^/(n) lj^.^(n) 


K‘^V{n) X 

where we used the Stirling approximation: 

lnU({A(n)}) ~A(n)/AlnA(n)/A . 
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In order to calculate the configurational partition function Z from Eq |13| we can introduce a saddle-point approxi¬ 
mation where the sum over all the distributions {A(n)} is replaced with the contribution from a distribution N*{n) 
which minimizes the free energy. We have: 


AT L ^ ^ V A‘‘*T’(n) A N ’ 


InZ 

~N~ 


(15) 
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and the distribution 7V*(n), that minimizes the system free energy, is: 


A^*(n) = No exp [ - ^Ae(n) - /iAu/(n)] . 


(16) 


where Nq is a normalization constant and ^ is a Lagrange multiplier associated to the constraint over the total free 
volume Vf. 

For simplicity we consider the case of hard spheres[22] where e(n) = 0. However the results obtained in this case 
can be extended to more general interactions, provided to substitute the free volume, with an effective temperature 
dependent free volume which takes into account the Boltzmann weight inside the cell|20j. To evaluate the sum in 
the continuum limit we must also introduce a function givf) such t hat q(vf)dvf gives the number of cells with free 


volume between Vf and Vf + dvf, assuming g{vf) ~ vj then from Eq, 
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the distribution of the cell free-volumes is 


p{vf) = 


N*ivM 


.,k-l 


N r(A:) (^yjY 


■ exp 




(17) 


where r(.) is the Gamma function, and k = ko + X with fco = d +1- Here A = 1 for low densities. Although the theory 
is valid for high densities the distribution Eqj^ works well also for small densities. For instance, it was shown in [23] 
that, in the zero density limit the Voronoi volumes for hard sphere in three dimensions follows Eq |17| with k = 5.586. 

In the case k = 1 the expression for p{vf) in Eq|17| reproduces the expected cell distribution of the free volume 
theory which is a pure exponential |5]. Moreover for fc> 1, Eq[l7| coincides with the so-called fc-Gamma distribution 
[14j that has been shown to reproduce well the distribution of the Voronoi volumes in granular packings for the entire 
range of densities investigated, with ko being a smooth function of the density. It was also pointed out in [ 13 ] that 


k = 



(18) 


where — (u/)^ is the variance of the free volume. 

This equation shows that in the free volume theory (fc = 1), where there is no cooperative diverging length , the 
fluctuation of the free volume vanishes as on the contrary if k diverges as fc ~ , the fluctuation of the 

free volume would vanishes as cr^ = (u/)^^^. Since, from Eqj^ k ~ and the compressibility kt it 

then follows kt -- ■ 

Later we will give arguments suggesting that ly = 1. Thus if we take the value of 9 = d/2 from RFTO, for d = 3, 
we have 


kt ~ ii/ro) , 


(19) 


relating the vanishing of the compressibility to the divergence of the cooperative length. 

Recently it has been shown that monodisperse disordered jammed particle are hyperuniform i.e. infinite wavelength 
volume fraction fluctuation vanish |24H27j with a structure factor S'(k) that tends to zero linearly in the wavenumber 
k as the glassy jammed state is approached, implying the vanishing of the compressibility kt S{0) and the 
divergence of a length ^dcf, relative to the direct pair correlation function p31 - EH] . More precisely since the Fourier 
transform of the direct pair correlation is given by c(k) = (<S'(k) — l)/pS'(k). Glose to the jamming glass state we have 
c(0) ^ —l/pS'(0), consequently, in 3d, since ^dcf = (~c(0))^'^^, we have ^dcf ~ ■ 

Interestingly from Eq[^ follows that the cooperative length ^ and the “hyperuniform” length ^dcf are related and, 
at the jamming glass state, ^ ^dcf- Note that our approach is only valid for equilibrium configurations, however 
it has been shown that for hard spheres, even at volume fraction not very close to jamming, the system appears to 
be out of equilibrium [55] 

Relaxation time. A density fluctuation can relax following two processes. The first is a local process related to 
the escape from the cage the second is a cooperative process, due to relaxation of a cooperative rearranging region. 

Escape probability. The escape probability associated with the first relaxation process can be computed by following 
the free volume theory m 


P — 

fisr. — 


poo 

/ P{vf) dvf , 

J V*r 


( 20 ) 


where p{vf) is given by Eq 16 and v*j is a free volume, below which the particle cannot escape from the cell. 
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The integral in Eqj^is: 


p — 

J- fifir. — 


nk,k^) 

m 


( 21 ) 


where r(., .) is the upper incomplete gamma function. 

Reconfiguration probability. The second mechanism that contributes to the relaxation time is the reconhguration 
of a cooperative rearranging region. The probability that such reconhguration occurs is given by Pcrr 1/^? where 
is the total number of conhgurations inside the droplet m, 0. In(flj) = where s 

- r 

Jo 


conhgurational entropy per particle, which from Eq 14 is Sc = B/X with B 


= In n is the 
’ p{vf)\np{vf) dvf being a slowing 


varying function of the density, which in the following we will consider constant with values of the order of 1. 
In conclusion, since A = we have for the conhgurational entropy 


9-d 


and for the reconhguration probability 


= Aexp(-B(^/ro)®) 


( 22 ) 


(23) 


where A is a constant. Eq|^is in agreement with the result of RFOT. 

The relaxation times associated with the two mechanisms are respectively proportional to the inverse of the prob¬ 
abilities: T ~ 1/Pesc and T ^ l/Pcrr- 

Relation between ^ and (vf). When the free volume per particle (vf) is small, even smaller than a particle can 
still manage to jump out of the cage albeit with an exponentially small probability, as the free volume distribution 
exhibits and exponential tail, which allows very rarely a single particle to have a free volume Vf > v*j. In order to 
have instead a cooperative relaxation of an entire droplet of radius R the free volume {R/rQ)‘^{vf) of the droplet must 
be large enough to allow the interior to change from one inherent state to another, keeping the boundary fixed. This 
transition can occur if the available free volume is distributed in such a way that all the free volume is essentially 
concentrated in a region around the surface of volume i?®, in order to match the internal inherent conhguration with 
the inherent configuration outside the boundary. The minimum size R — ^ must be such that the particles in the region 
near the boundary has just the minimum free volume, roughly of the order of v*, enough to realize the matching. This 
condition corresponds to the marginal stability limit [16]; namely, for fixed boundary, droplets of radius R < ^ are 
stable, while droplets with radius R> ^ are unstable. In conclusion this condition implies (C/?’o)‘^(u/) = ki(^/ro)^Vj, 
where ki is a constant, from Eq{n] we obtain: 


A = kiVf/{vf). 


(24) 


We expect ki < 1 because A must become of the order of 1 when the average free volume (vf) is larger than Vj. 

Using the value 9 = d/2 and (^/ro)‘^/^ = kiv‘f/{vf) it follows k = ko + kiv*f/{vf). Note that (Eq 24) is valid only 
for A > 1. 

To calculate the relaxation time we have to distinguish two regimes: i) a low density regime where A = 1; ii) an 
high density regime where A > 1. 

1) In the low density regime, where A = 1, the relaxation time is only given by the escape process. From Eq |21| we 
have: 


tq 

p 

J p.sf 


(25) 


with Pesc given by Eq |21| with k = ko + 1. 

The incomplete Gamma function can be evaluated for small and large densities. For low densities, when 
(vf) > Vp Pesc decreases and apparently goes to zero following a power law behaviour: 


To 


{{vf)/v}-b) 


-7 


(26) 


with b close to 1 and 7 = 1 . This apparent power law is only a crossover towards an exponential behavior for 
high densities {vf) < vy. 


To 


exp {{ko + l)vy{vf)) 


(27) 
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FIG. 1: Demonstration that cell theory predicts well experimental measurements of the relaxation time. The figure reports 
the logarithm of the relaxation times r vs. relative inverse mean square rattling amplitude (u?) / (w^) for a large number of 

The blue dotted line and the 
contributes significantly only 


different glass-former materials from (Fig.3). The green line is the behaviour given by Eq 28 
green dotted line are respectively the contribution from and Pcrr- One can note that P^ 
at very high densities ((ug) / > 1). 


2) In the density regime, where A > 1, both the escape and collective rearrangement mechanisms contribute to the 
relaxation time. The relaxation time is given by 


— — l/(Pesc + Pcrr) 
To 


(28) 


with Pgsc given by Eqj21l Pcrr given by Eq 23 and k = ko + kiv*j:/{vf). This gives an expression for the relaxation 
time as function of VjJJvf). 

Using the expression for the incomplete Gamma function evaluated at large values of v*^/{vf) the relaxation time at 
high densities is given more explicitly by 


■^0 exp (^-kov*j:/{vf) - ki{vy{vf)Y'^ + Aexp (^-kiBvy{vf) 


(29) 


In a recent work by Larini et al. |28j it was pointed out that a large numbers of glass-former materials and several 
polymers reveal a universal correlation between the structural relaxation time (and the viscosity) and the mean square 
rattling amplitude (or Debye-Waller factor), of the vibrational modes. It is reasonable to conjecture that the 

mean square rattling amplitude must be related with the free volume, with larger amplitudes associated to greater 
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free volume availability. On the basis of dimensional analysis, we expect (u^) / (ug) = {{vf) with and 

Vg respectively the mean square rattling amplitude and the average free volume at glass transition. By using the 
experimental data collected by Larini et al. in [2H] (Fig-3) we can compare the cell theory expression for r from Eqj^ 
with the experimental data. The result is shown in Figj^were we can see that Eq |28| reproduces well the experimental 
findings with a htting curve (green line) with parameters: fco = 0.68, ki = 0.16, B = 8.24, tq = 0.86, A = 0.018 and 
Vg = 0.12v}. 


Finally note that if we use, from the free volume theory, the relation between free volume and temperature (vf) 

T — Tq and 9 = d/2 we recover the RFOT expression for ^ ~ (T — (Eq 241, Sc ~ (F — Tg) (Eq 221 and, from 

Eq |29[ a VET for the relaxation time at a temperature very close to the Kauzmann temperature Tg. 

In conclusions using a cell theory previously developed, combined with the Adam-Gibbs-Wolyness cooperative 
rearranging region approach, we have shown that there are two distinct regimes characterizing the relaxation time. 


In the first regime, where the cooperative length ^/rg = 1, the relaxation time is entirely due a local process Eq 25 


This regime is characterised by an apparent power law, reminiscent of the mode coupling theory behaviour, followed 
by an exponential behaviour, typical of the free volume theory. At higher densities the relaxation is due to both escape 
process (Eq{2T|) and relaxation of cooperative regions (Eq|23|. These two processes coexist until at high density the 
presence of the cooperative regions suppress the escape probability and the only process left is due to the rearranging 
cooperative regions. For the relaxation due to the local process we use the relation which comes from the free volume 
theory. This is a simplification as the local relaxation may be more complex due not only to the free volume but to 
other quantities, involving also second and third neighbors!^. 

Interestingly our approach relates also the cooperative length ^ and the hyperuniform uniform length ^dcf ■ 
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